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ABSTRACT 

We have developed a physically self-consistent model of the disk around the 
nearby 10 Myr old star TW Hya which matches the observed spectral energy 
distribution and 7mm images of the disk. The model requires both significant 
dust size evolution and a partially-evacuated inner disk region, as predicted 
by theories of planet formation. The outer disk, which extends to at least 140 
AU in radius, is very optically thick at infrared wavelengths and quite massive 
(~ O.O6M ) for the relatively advanced age of this T Tauri star. This implies 
long viscous and dust evolution timescales, although dust must have grown 
to sizes of order ~ 1 cm to explain the sub-mm and mm spectral slopes. In 
contrast, the negligible near-infrared excess emission of this system requires that 
the disk be optically thin inside <; 4 AU. This inner region cannot be completely 
evacuated; we need ~ 0.5 lunar mass of ~ 1 /iin particles remaining to produce 
the observed 10/im silicate emission. Our model requires a distinct transition in 
disk properties at ~ 4 AU, separating the inner and outer disk. The inner edge 
of the optically-thick outer disk must be heated almost frontally by the star to 
account for the excess flux at mid-infrared wavelengths. We speculate that this 
truncation of the outer disk may be the signpost of a developing gap due to the 
effects of a growing protoplanet; the gap is still presumably evolving because 
material still resides in it, as indicated by the silicate emission, the molecular 
hydrogen emission, and by the continued accretion onto the central star (albeit 
at a much lower rate than typical of younger T Tauri stars). TW Hya thus may 
become the Rosetta stone for our understanding of the evolution and dissipation 
of protoplanetary disks. 
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1. Introduction 

The discovery of extrasolar planets (Marcy & Butler 1998 and references therein) has 
opened up a new era in the study of planetary systems. While many important clues to the 
processes of planet formation can be obtained from studies of older systems, the best tests 
of formation scenarios will require the direct detection of actively planet-forming systems. 

It is thought that the formation of giant planets involves the sweeping up of material 
in a wide annulus in the circumstellar disk, resulting in the development of a gap (Lin 
& Papaloizou 1986, 1993; Bryden et al. 1999). Material inside the planet-driven gap can 
continue to accrete onto the central star; if the planet can prevent material from accreting 
across the gap into the inner disk, the eventual result would be the evacuation of the region 
interior to the planet. In the case of the Solar System, the formation of Jupiter might have 
prevented outer disk gas from reaching the inner solar system; the inner gas disk accreted 
into the Sun, while solid planetesimals remaining behind eventually formed the terrestrial 
planets. 

The above scenario suggests that the signature of a forming giant planet would be 
the presence of a gap that is not entirely evacuated. In this case dusty emission from the 
inner disk might still be observable. In addition, giant planet formation requires the prior 
consolidation of large solid bodies to serve as cores for subsequent gas accretion; while 
such bodies would be invisible with current techniques, one would expect to see evidence 
for substantial growth in dust particles. Finally, if the inner disk has not been completely 
evacuated by the forming planet, one might expect to observe continued accretion onto the 
central star, as all T Tauri systems with inner disks as detected from near-infrared disk 
emission are also accreting (Hartigan et al. 1990). 

In this article we propose that the relatively young low-mass star TW Hya has a 
developing gap in its inner disk qualitatively similar to that expected from planet formation. 
The evidence supporting this this proposal is: (1) reduced (optically-thin) emission from 
the inner disk; (2) mm-wave spectra which seem to require grain growth; (3) extra emission 
from the edge of the outer disk; and (4) continued accretion onto the central star, albeit at 
a rate substantially lower than that observed from most T Tauri stars. TW Hya has an age 
of 10 Myr and so according to current theories it is quite likely to be close to the epoch 
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of planet formation. It is also part of an association of young stars of similar age, but it 
stands out in that it is the only system still accreting at a substantial rate (Muzerolle et al. 
2000). The (outer) disk of TW Hya is quite massive, so that there is likely to be more than 
enough material available to form giant planet (s). TW Hya is also the closest known such 
system, and so it will be a prime target for following studies to confirm our model. 

2. Observations 

The observations used to constrain our disk model are taken from the literature. In 
addition, we have obtained a narrow-band L (3.55-3.63 microns) measurement of TW Hya 
on June 14th, 2000 with the Australian National University 2.3-m telescope at Siding 
Spring, using the near infrared camera CASPIR (Cryogenic Array Spectrometer Imager; 
McGregor et al. 1994). CASPIR contains a 256x256 InSb array. TW Hya was imaged at L 
with a pixel scale of 0.25 arcseconds and an on-source integration time of 236 seconds. The 
standard star BS4638 (magnitude 4.50 at L) was observed at a similar airmass to TW Hya 
and was used to calibrate the magnitude of TW Hya. The L magnitude we obtained was 
7.12, equivalent to 0.39Jy. 

This measurement is important in constraining the amount of near-infrared excess in 
the system. Using K=7.37 (Webb et al. 1999), K-L = 0.25. The K-L color for a K7V star 
is 0.11 (Kenyon & Hartmann 1995), which would yield an infrared excess of 0.14. However, 
such small infrared excess does not necessarily imply emission from disk material (Wolk 
& Walter 1996). For example, the non-accreting stars in Taurus (spectral types ~ K7 to 
M2) have K-L between -0.05 and 0.25, while 95% of the accreting stars in Taurus have 
K-L between 0.35 and 1.2 (Meyer, Calvet, & Hillenbrand 1997). This small value of K-L is 
consistent with the simultaneous flux measurements of Sitko et al. (2000) indicating little if 
any hot dust emission at wavelengths < 5fim, though it is difficult to rule out completely 
any excess. 

3. Model assumptions 

We assume that the heating of the disk is due to stellar irradiation and viscous 
dissipation, and calculate self-consistently the disk heights and temperatures following the 
methods we have applied to interpret young ~ 1 Myr old disks (D'Alessio et al. 1998 = 
Paper I, 1999 = Paper II, 2000 = Paper III.) We adopt M„ = 0.6 M , R* = 1 R and 
T* = 4000 K for the stellar mass, radius, and effective temperature (Webb et al. 1999), 
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and a mass accretion rate of M = 5 x 10~ 10 M yr~ 1 as derived by Muzerolle et al. (2000), 
which we take as constant through the whole disk. The inclination of the disk axis to the 
line of sight is i ~ 0°, in agreement with the nearly symmetric HST images (Krist et al. 
2000), and we adopt the Hipparcos distance of 55 pc (Wichmann et al. 1998). 

With the assumption of steady accretion, the disk surface density scales as 
E oc Ma' 1 !" 1 , with T is the midplane temperature (see Paper II). The temperature 
does not vary much between acceptable models, partly because it is mostly determined 
by irradiation heating, not viscous energy dissipation. The disk mass therefore is roughly 
constant for constant Ma' 1 . Moreover, for long wavelengths where much of the disk is 
optically thin, the emergent fluxes also tend to scale in the same way. For the purposes of 
fitting, we vary the parameter a as the dust opacities are varied, but note that varying a 
is equivalent to varying the disk mass, because M is fixed. To first order, however, choices 
with Ma" 1 ~ constant are acceptable. 

We use a dust mixture consisting of silicates, refractory organics, troilite and water ice, 
following Pollack et al. (1994 = P94). The dust size distribution is taken to be n(a) ~ a~ p , 
with p = 3.5, between given a m i n and a max . Optical properties for the compounds are taken 
from Jager et al. (1994), P94, Begemann et al. (1994; also see Henning et al. 1999), and 
Warren (1984). We consider the grains to be compact segregated spheres, and calculate the 
opacity using a Mie scattering code (Wiscombe 1979). Sublimation temperatures for the 
different grain types are taken from P94. 

4. Disk model 

The spectral energy distribution (SED) of TW Hya is shown in Figure 1. References 
for the observational data are given in the figure caption. Note that our narrow-band L 
observation overlaps with Sitko et al. (2000) fluxes. 

In Figure 1, we compare the SED with the median SED for ~ 1 Myr-old T Tauri stars 
in the Taurus molecular cloud (Paper II), normalized to the TW Hya stellar photosphere 
at H (1.6yum), thus compensating for the differing distances and stellar luminosities. Even 
though the fluxes of TW Hya are relatively high compared with the median Taurus SED 
at A ^ 20/iin, there is a large flux deficit in TW Hya below 10 microns (Jayawardhana et 
al. 1999); in particular, fluxes are essentially photospheric below ~ 6/im (Sitko et al. 2000). 
The flux deficit below 10 /xm has lead to inferences of disk clearing inside a few AU from 
the central star (Jayawardhana et al. 1999), but since the disk is still accreting mass onto 
the star (Muzerolle et al. 2000), it has to extend all the way into the corotation radius at 
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least, so the inner disk radius has to be ^ 0.03 AU (inferred from the 2 day photometric 
period [Mekkaden 1998]). 

We show below that models with a uniform dust well-mixed with the gas throughout 
the disk extending all the way onto the inner radius cannot explain these features of the 
SED. The observations are much better understood if the disk is divided into two regions: 
the outer disk, which is more nearly comparable to the structure inferred for typical T Tauri 
disks (Paper III), and the inner disk, which is much more optically thin than in typical T 
Tauri disk models. 



4.1. Outer disk 

As discussed in Paper II, ISM dust mixtures (with small a max ) cannot explain the 
far-IR and mm-wave fluxes of T Tauri stars. The similarity between the median Taurus 
SED and that of TW Hya at wavelengths A ^ 100 /im suggests disk models of the type 
explored in Paper III, in which we allow for growth to large particles, can in principle 
explain the observations. 

Figure 2 shows results for flared, irradiated disk models calculated with the methods 
of Paper III. The model SEDs shown are for a max = 1 nim (which fits the median SED 
in Taurus; Paper III), 1 cm, and 10 cm, calculated with abundances in the dust mixture 
usually assumed for protoplanetary disks (P94), which yield a dust-to-gas ratio of 0.013. 
It can be seen that models where grains have grown to a max ~ 1 — 10 cm provide a much 
better fit to the long-wavelength SED than the a max = 1 mm model, requiring no or very 
little additional emission from a wind or non-thermal sources. Models with a max <C 1 mm 
fail to reproduce either the sub-mm and mm spectral slopes or the total flux levels for 
reasonable disk masses. 

Assuming an outer disk radius of ~ 140 AU, comparable to the radius at which Krist 
et al. (2000) find a rapid decline in disk density, the disk mass is ~ 0.03, 0.06, and 0.1 1M 
for the a max = 1 mm, 1 cm, and 10 cm models, or a = 5 x 10~ 4 , 3 x 10~ 4 , and 1 x 10~ 3 , 
respectively. The a max = 1 cm model has a Toomre parameter Q ~ 1 at its outer edge and 
so is near the limit expected for gravitational stability. The high mass values are due to the 
opacity at 7 mm, which decreases as a~^J x 2 (for p = 3.5, see Paper III), so higher masses 
are needed to account for the flux for larger grain mixtures. For example, for the a max = 
10 cm mixture, the opacity at 7 mm is a factor of 3 lower than the frequently assumed law 
k>bs — 0.1(A/250 fxm)^ 1 (Beckwith & Sargent 1991) and the slope is slightly flatter, oc A - 8 
(cf. Paper III, Figure 2). 
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Our total disk mass estimates, using the dust-to-gas ratio of 0.013, are much higher 
than the gas mass obtained by Kastner et al. (1997) from 12 CO emission, 3.5 x 10~ 5 M Q . 
This discrepancy may be attributed to a large degree of molecular depletion or to the 
fact that optical depth effects may not have been properly included (Beckwith & Sargent 
1993); it could also be due to molecules existing in the gas phase only in the hot upper 
atmospheric layers of the disk where only a small amount of mass resides (Willacy & 
Langer 2000). An alternative to high disk masses, not considered in this work, is to have 
larger opacities at 7 mm than we obtain; porous aggregates, specially for amorphous carbon 
particles, may result in such larger opacities (Stognienko, Henning & Ossenkopf 1995). Our 
results for grain growth are not unique. However, within the assumptions we have made 
concerning dust opacities, we cannot reproduce the mm-wave fluxes and spectral slope 
without including dust particles much larger than those of a "standard ISM" mixture (see 
Paper III). Beckwith & Sargent (1991) noted that optical depth effects could make the 
mm-wave spectral slope flatter, and thus reduce or eliminate the need for grain growth. 
However, we are unable to make TW Hya disk sufficiently optically thick over its large 
radial extent. 



4.2. Edge of the outer disk 

The well-mixed grain-growth models of Figure 2, with self-consistently calculated 
temperature structures, exhibit too little emission in the 20 — 60/im wavelength region in 
comparison with observations. In principle, disk models with more small particles could 
have higher fluxes. However, it is apparent that these models also predict far too much 
flux at wavelengths ^ 10/im, because they are too optically thick in the inner regions; this 
suggest that some clearing has occurred in the inner disk. Moreover, it suggest that the 
outer disk should have an optically-thick edge in which extra heating could be important. 

In the case of irradiation of an optically thick disk with a smoothly-varying thickness, 
the stellar flux captured by the disk (and thus the disk heating per unit area) depends 
on the cosine of the angle between the direction of incidence of the stellar beam and the 
normal to the disk surface /io (Kenyon & Hartmann 1987). In general, stellar radiation 
penetrates the disk very obliquely, so /io is fairly small, ~ 10 -3 — 10 -2 (cf. D'Alessio 1996). 
However, if the disk had an inner edge, this portion of the disk would be illuminated by the 
star more directly, increasing /io dramatically and thus increasing the amount of irradiation 
heating. We propose that the outer regions of the TW Hya disk can be described as in §4.1, 
but that it is truncated at a few AU by a steep, optically thick region, where most of the 
mid-IR flux excess arises (Figure 3). Inside this region lies the inner optically thin disk, 
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which produces negligible continuum flux. 

To calculate the emission of the disk edge, we have assumed that its irradiation surface, 
that is the surface where most of the stellar energy is deposited, has a shape given by 
z s = z (R )exp[(R — R )/ AR], where z Q is the height of irradiation surface of the outer disk 
at radius R Q and AR is a characteristic width of the edge. The edge continuum emission 
is produced in the photosphere, which for simplicity we take as the irradiation surface; 
above this layer, there is a hotter optically thin region, the atmosphere, which we take as 
isothermal (cf. Chiang & Goldreich 1997). We think that a more detailed treatment of the 
structure of this region is not necessary given the exploratory nature of this study. The 
temperature of the edge photosphere, assuming that half of the intercepted flux reaches the 
photosphere, is given by 



where /xo = cos(9 ) is obtained from O = tt/2 — tan 1 (dz / dR). The flux is evaluated as 

F v = f R ° I u 27cRdR, (2) 

JRth 

where R t h is a the radius where the disk becomes optically thin, and I v = B u {T phot ). 

The location of the outer disk edge and its width, characterized by R Q and ~ AR, are 
fairly well restricted by the mid-infrared excess, although the actual shape is not so well 
constrained as long as most of it faces the star. In general we find AR ~ 0.5 AU, because 
the shape of the excess is fairly narrow and thus cannot be produced by a region with a 
large range of temperatures. In addition, if R Q is much smaller than ~4 AU, then the edge 
is too hot and there is too much excess below 10 fim. Similarly, if R Q is much larger, there 
is too little flux. With R Q ~ 4 AU and AR ~ 0.5 AU, /i G varies from ~ 0.2 at R = 3 AU 
to ~ 0.7 at R = 4 AU, and T p f lot varies from ~ 80 to 100 in this range. The resulting 
continuum flux is shown in detail in the lower panel of Figure 4, while the upper panel 
shows the fit to the SED of the composite disk model. 

The temperature of the optically thin atmosphere of the edge is given by 



where W(R) is the geometrical dilution factor, W(R) = Q^/Atc « (R*/2R) 2 , for R » R*, 
and k* p and Kp are Planck mean opacities calculated at the stellar and local radiation fields, 
respectively (cf. Paper I, Paper II, Paper II). The contribution to the flux from this region 
is given by eq. (|), with I v = B u (T up )t u and fx*-, where r* and x* are the optical 

depth and extinction coefficient at the characteristic wavelength of the stellar radiation. We 




(1) 



W{R)n P Tt = K P {T up )T* v 



(3) 
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follow Natta et al. (2000) taking r* pa /i , but in this work we include the effect of scattering 
at the wavelengths at which the stellar radiation is absorbed. The neglect of scattering in 
the calculation of x* leads to artificially large values of r u and thus of I u , and would be 
only appropriate if scattering was completely forward. However, the asymmetry parameter 
g = (cosO), with the scattering angle, is in general < 1, and this approximation is not 
valid. We have included the effect of scattering using x* = K * + (1 — g*)cr*, where fi^*, c*, are 
the absorption and scattering coefficients and the asymmetry parameter for the assumed 
dust mixture, all evaluated at the characteristic wavelength of the stellar radiation. 

We assume that the atmosphere of the edge has small grains so it can produce emission 
in the silicate feature. We considered different glassy and crystalline pyroxenes and olivines 
(with optical properties from Laor & Draine 1993, Jager et al. 1994, and Dorschner et al. 
1995). We find that glassy pyroxene Mg - 5 Fe 0A3 Ca°- m Al 0m SiO 3 has one of the highest 
ratios as required for producing a strong band (see also Natta et al. 2000). However, 

even with the increasing heating at the disk edge, the conspicuous silicate emission feature 
seen in the SED of TW Hya cannot be explained by emission from the atmospheric layers of 
the edge. This is due in part to the inclusion of scattering of stellar light in the calculation 
of x*i as discussed above, and in part to the adopted dust mixture; organics dominate the 
opacity at the stellar wavelength (Paper III, Figure 2), yielding Kio^m/x* ~ 1- Note that 
the atmosphere of the optically thick outer disk, which is irradiated less frontally by the 
star and is thus cooler than the atmosphere of the edge, cannot produce significant silicate 
emission either. 



4.3. Inner disk 

Since the region inside the outer disk edge is not empty because material is still 
accreting onto the star (Muzerolle et al. 2000), we have explored the possibility that a small 
amount of particles coexists with the accreting gas of the inner disk, giving rise to silicate 
feature emission, while on the other hand still resulting in negligible continuum flux. 

We assume that the temperature of the dust in this region is given by equation (|3|). The 
emergent intensity is calculated as I u = B u k u t w / ^(lOyum), where r 10 , the optical depth at 
10/im, is a parameter. This region extends to ~ 0.02 AU, the radius where grains sublimate. 
We find that we can fit the silicate feature with t\q ~ 0.05, see Figure 4, lower panel; this 
optical depth is small enough for the disk continuum emission in the near infrared to be 
negligible. The best fit to the profile is achieved with glassy pyroxene, in agreement with 
the fact that the profile is closer to that in young stars than in comets (Sitko et al. 2000). 
The grain sizes are in the range a m i n ~ 0.9/im and a max ~ 2/zm. Smaller grains have much 
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higher temperatures and produce too much emission at short wavelengths, resulting in a 
narrower silicate profile than observed. Bigger grains produce too little emission. 

From Tio ~ 0.05, we get a column density for the ~ l/im dust of S d ~ 4 x 10 _3 gcm~ 2 , 
which implies a mass inside 4 AU of ~ 4 x 10 25 g ~ 0.5 lunar masses. A strong lower 
limit to the mass of gas in the inner disk can be obtained from the observed rate of gas 
accretion onto the star assuming that the material in the inner disk is in free-fall towards 
the star, M gas » MR 3 / 2 /(2GM*) 1 / 2 . Inside 4 AU, we obtain M gas » 6 x 10~ 10 M Q . 
Another lower limit can be obtained from S^, assuming a normal dust-to-gas ratio, yielding 
S g ~ 0.4gcm~ 2 and a mass M gas > 2 x 1O _6 M ~ 0.6 earth masses. If solid material is 
hidden into large bodies, then the mass in gas could be much higher than this value. From 
this limit, nonetheless, we can get an upper limit for the radial velocity of the gas at 1 
AU, v R < MjlixKLg ~ 0.008 km s" 1 << v K (lAU) ~ 23kms _1 , so the gas probably drifts 
inwards slowly, following nearly Keplerian orbits in the inner disk. 

Weintraub, Kastner, and Bary (2000) have detected a flux of ~ 1.0 x 10~ 15 ergs _1 cm -2 
in the 1-0 S(l) line of H 2 towards TW Hya. Using Tine et al. (1997) models, we can estimate 
a flux in this line at 55 pc, F(H 2 ) ~2x lO- 15 (e/lO- 2 V 1 )(M H2 /lO- 8 M )ergs~ 1 cm" 2 , 
where e is the emissivity per H 2 molecule. For a gas temperature of ~ 1000K and a 
ionization rate of 3 x 10 _10 s _1 , estimated with typical parameters and an X-ray luminosity 
L x ~ 10 30 ergcm~ 3 (Kastner et al. 1997) using the Glassgold, Najita, & Igea (1997) 
corrected expression, e is 2.3 x 10 _22 ergs _1 for densities nn > 10 7 cm~ 3 (because the line 
becomes thermalized; S. Lepp, personal communication). Assuming a disk height of 0.05 
AU, we estimate tin » 4 x 10 7 cm~ 3 inside 4 AU, so from the mass limit estimated above, 
we see that the observed H 2 very likely arises in the inner disk. Firmer conclusions require 
a determination of the gas temperature in the inner disk, which is left for future work. 



4.4. Model tests: comparison with VLA 7 mm Data 

In principle, observations of disk surface brightness distributions can be used to test 
our physically self-consistent models based on SED fitting. The optical and near-infrared 
scattered light distributions presented by Krist et al. (2000) and Trilling et al. (2001) do 
not have sufficient resolution to probe our inner disk region; however, they do provide 
constraints on our outer disk structure. We find that the scattered light fluxes predicted 
by the a max = 1 cm model presented above matches the observations of Krist et al. (2000) 
reasonably well, although Krist et al. find evidence for structure that cannot be reproduced 
in detail with disk properties smoothly-varying with radius, as assumed here. This work 
will be reported in a forthcoming paper (D'Alessio et al. 2001). 
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The high angular resolution data at 7 mm from the Very Large Array (Wilner et al. 
2000) begin to approach the resolution needed to probe the inner disk, as well as constrain 
the thermal dust emission of the outer disk instead of simply the scattering surface, and 
thus provide an additional test of our model. Moreover, the imaging data provide an 
independent test, since the model was constructed only to match the SED. The value of no 
sets an upper limit to the opacity at A > 10/xm, so the inner disk should not contribute 
significantly to the emission at 7 mm. We have calculated the intensity profile of the disk, 
assuming different values for the brightness temperature of the edge between 3 and 4 
AU, since we do not know the properties of its interior. An upper limit for is the surface 
temperature, ~ 100K, but it could be lower if it was optically thin at 7 mm, so we have 
varied T b between 60K and 100K. 

Figure [5]a shows the 7 mm images at two angular resolutions. The "low" resolution 
image (~ 0'.'6, about 35 AU) emphasizes extended low brightness emission, and the "high" 
resolution (~ O'/l, about 6 AU) shows the smallest size scales where emission remains 
detectable with the available sensitivity. Figure |5|b shows the result of imaging the model 
brightness distribution at these two angular resolutions using the same visibility sampling as 
the observations. For a more quantitative comparison, Figure [|c shows the residual images 
obtained by subtracting the model visibilities from the 7 mm data, and then imaging the 
difference with the standard algorithms. 

The low resolution model image, with the centrally peaked dust emission (T& > 60K), 
compares very favorably with the 7 mm images of Figure |^a; unfortunately, the higher 
resolution is still not adequate to resolve directly the presence of the inner hole. But in 
both cases, the residual images show only noise. The agreement between the model and 
observations confirms the insightful parametric modeling of Wilner et al. (2000), who fitted 
the observations with radial power laws for the temperature and surface density. This 
was to be expected, since Wilner et al. adopted powers that were consistent with the 
predictions of irradiated accretion disk models at large radii. However, since our models are 
constructed from first principles, by solving the disk structure equations in the presence of 
stellar irradiation subject to the constraint of steady accretion, we are able to predict and 
confirm not only the radial dependences of the physical quantities (which are not exactly 
power laws), but also their absolute values. The free parameters of our model are essentially 
the disk mass (or equivalently, the parameter a, see §3) and the dust mixture. The models 
were constructed primarily to fit the SED, but they can also fit the radial brightness 
distribution, as the good agreement between the observations and the model indicate. On 
the other hand, although there is no guarantee that the outer disk is completely steady, or 
has a constant a, or ours is the correct prescription for the dust mixture, these assumptions 
seem to provide a good description of the physical situation in the disk of TW Hya. 



5. Discussion 



The lack of infrared flux excess at wavelengths below 10 microns could be easily 
understood if the inner disk was evacuated inside a few AU. However, the fact that material 
is still being transfered onto the star, as evidenced by the broad emission line profiles and 
the ultraviolet excess emission (Muzerolle et al. 2000), indicates that gas still exists in 
the inner disk. However, this material must be optically thin to explain the flux deficit 
in the near infrared. If the outer disk extended inwards to the magnetospheric radius, 
then the surface density S at 2 AU would be ~ 320gcm" 2 , implying a surface density 
in solids of ~ 3gcm -2 (for a dust to gas mass ratio of 0.01). To obtain an optical 
depth at the near-infrared of < 0.05 with this amount of dust, the opacity should be of 
order Kdust < 0.05/Sd ust ~ 0.02cm 2 g _1 . If we assume that the grains are large enough that 
the opacity per gram of grain is ~ Qext^c 2 / (4/3)7rp ff a 3 , with Q ext the extinction cross 
section and p g the grain density, both or order 1, then K dust ~ 1/a. Using this estimate, we 
find that the original solids should have grown to sizes a > 50 cm to account for the low 
optical depth in the near infrared. While this estimate is very rough, it suggests that very 
significant grain growth has probably already taken place in the inner disk of TW Hya, 
especially if the solid to gas ratio has increased significantly from solid matter left behind 
as gas flows onto the star. This inside-out coagulation and settling is consistent with the 
evolution of solid particles predicted in solar nebula theories. Dust grains are expected to 
coagulate by sticking collisions at a rate oc 1/Q 2 oc P 2 and reach a maximum size at the 
midplane in a few xlOOOP, where Q is the Keplerian angular velocity and P the orbital 
period (Safronov 1972; Goldreich & Ward 1973; Weidenschilling 1997; Nakagawa 1981). 
Since P ~ 1.4 x 10 3 (i?/100A[/) 3 / 2 yrs, coagulation and settling are expected to occur first 
in the inner disk. 

In the outer disk, a dust mixture where grains have reached ~ 1 cm and dust and gas 
are well mixed can explain the observations, within the framework of our assumptions. This 
interpretation is consistent with the Weidenschilling' s (2000) argument that grain growth 
to at least 1 cm is necessary before particles start settling to the midplane, because smaller 
grains get stirred around by turbulence, and we suggest that turbulence plays an important 
role in the evolution of dust in protoplanetary disks. Grains in the upper layers can absorb 
stellar radiation and efficiently heat the disk and make it flare, resulting in the large far 
infrared and millimeter fluxes observed. Nonetheless, grains in the outer disk have grown 
significantly from ~ 1 mm, the size that characterizes the median SED of Taurus (Paper 
III). To get high mm fluxes for such large grains, the disk mass has to be high, ~ O.O6M , 
since the opacity at 7 mm decreases as grains grow (cf. Paper III, §2). This high value for 
the mass is surprising; it is similar or higher than values typically found for disks around 1 
Myr old T Tauri stars (Beckwith et al. 1990), and comparable to that of objects emerging 
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from the embedded phase (D'Alessio, Calvet, & Hartmann 1997). On the other hand, as we 
discussed in Paper III, disk masses may have been underestimated if the standard power 
law opacity is applied indiscriminately to disks where significant grain growth may have 
taken place. 

It may be possible that the 3.6 cm data point contains a significant contribution from 
hot (ionized) stellar plasma. This contribution is unlikely to come from a wind. If the 
flux emitted by the wind is proportional to the mass loss rate, then scaling down from 
the wind/jet contribution in the case of HL Tau, which has M ~ lO^M^r- 1 to the 
case of TW Hya, and assuming a ratio of mass loss rate to mass accretion rate of ~ 0.1 
(Calvet 1998), we find that this contribution at 3.6 cm should be ~ 2 x 1CT 4 mJy, too 
low to be important. A non-thermal contribution due to stellar surface activity is still 
possible, given the high X-ray flux of TW Hya (Kastner et al. 1997). Second epoch 
observations should resolve this question, because any non-thermal contribution should be 
highly time-dependent. In any event, the high flux at 7 mm, at which a non-thermal plasma 
contribution is much less likely, indicates that the grain mixture must contain grains that 
have grown at least to ~ 1 cm, otherwise the opacity would already turn to the optical 
limit and drop, as is apparent in Figure 2 (cf. Paper III, Figure 2). 

The peculiarities of the disk of TW Hya suggest an advanced state of evolution of 
the dust. In addition, the presence of a steep transition between the inner and outer disk 
suggest the action of an agent other than dust settling to the midplane, which in first 
approximation would be expected to produce effects that vary monotonically with radius. 
We speculate that we may be seeing the outer edge of a gap opened by the tidal action of a 
growing protoplanet. Numerical models of gap formation in disks result in a sharp drop of 
the surface density at the gap and enhancements near the edges, as disk material is pulled 
away from the protoplanet and into the outer disk (Bryden et al. 2000; Nelson et al. 2000). 
The edge of the gap would be facing the star, and if it was optically thick, would result 
in excess emission similar to that modeled schematically by our truncation region. The 
situation could be comparable to the one considered by Syer & Clarke (1995), who studied 
the evolution of the disk spectrum under the presence of a planet formed in the disk, in 
the case in which the mass of the protoplanet was larger than the mass of the disk at the 
planet location. In this case, the disk cannot push the planet inwards at the local viscous 
time-scale and material forms a reservoir upstream of the gap (Lin & Papaloizou 1986). 
The surface density enhancement at the edge could make it optically thick, so the midplane 
temperature would not increase as much as the surface temperature; the scale height, in 
this case, would not increase significantly, justifying our assumption of similar heights for 
the edge region and outer disk and the neglect of the effects of shadowing. The details of 
the structure of this region, in any event, are of necessity very schematic. The extrapolated 
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mass of the disk inside ~ 4 AU is ~ 1(T 3 M ~ 1 Jupiter mass, so this mechanism would 
work if a giant planet were forming in this region of the disk. 

Finally, we address the question of what TW Hya tells us about protoplanetary disk 
evolution. Although we do not have direct measurements of the gas mass of the outer 
disk, it would be very surprising if the gas were strongly depleted there, as suggested by 
Zuckerman et al. (1995). Both the SED and the scattered light emission as a function 
of radius follow the predictions of models in which the dust is well-suspended in the gas, 
producing finite scale heights of the dust disk. Based on this inference, we suggest that the 
existence of the (fairly massive) TW Hya disk is evidence that neither UV radiation nor 
stellar winds can efficiently remove gas at radii ^ 5 AU over timescales of 10 Myr for stars 
formed in dispersed, isolated environments. 

The survival of disk material at 10 Myr also places constraints on viscous timescales. 
The surface density of an accretion disk evolves with time, as viscous dissipation makes the 
disk transfer mass onto the star and expand to conserve angular momentum in a viscous 
time-scale given by -R 2 /V, where R c is a characteristic radius and v the viscosity (cf. Pringle 
1981). Here we use the similarity solutions for disk evolution of Hartmann et al. (1998), 
which assume temperature and surface density distributions comparable to our outer disk 
models. The viscous time-scale in these solutions is t s ~ R\jv ~ 8 x 10 4 (10~ 2 /a)(-Ri/10AC/) 
yrs, where R\ is the radius that contains ~ 60 % of the initial mass (at t = 0), and the 
viscosity v is expressed in terms of the standard Shakura-Sunyaev (1973) a parameterization. 

Adopting q~5x 10~ 4 and R\ ~ 100 AU, then i s ~ 1.6 x 10 7 yrs. In the similarity 
model, the mass of the disk falls in time as Md oc T -1 / 2 , where T = 1 + t/ts. Thus, 
if the total disk mass at 10 Myr is 0.06M Q , then the initial disk mass must have been 
0.08M Q , comparable to that of objects still surrounded by infalling envelopes, like HL 
Tau (D'Alessio, Calvet, & Hartmann 1997). The surface density drops below a S oc R~ l 
power-law by a factor of 1/e at a distance Ri(T) = RiT ~ 160 AU, close to the ~ 140 AU 
transition in midplane density found in Krist et al. (2000). With these values of disk mass 
and radius, the value of the Toomre parameter Q remains close to unity and so gravitational 
instabilities (which could transfer angular momentum rapidly and thus invalidate the purely 
viscous evolution model) need not arise. 

The mass accretion rate for this constant-a model Q at an age of 1 Myr is 
M = Md(0) I (2t s T 3//2 ) ~2x 10~ 9 MQyr -1 , near the lowest values observed in accreting T 
Tauri stars of this age (Gullbring et al. 1998; Hartmann et al. 1998). The mass accretion 
rate at 10 Myr is then ~ 1 x 10~ 9 Mqi/t^ 1 , which is slightly higher than estimated by 



6 Note that a factor (Md(0)/0.1M Q ) is missing in expression (39) of Hartmann et al. (1998). 
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Muzerolle et al. (2000). It would not be surprising if accretion is being halted in the inner 
disk, and that inner and outer disk are becoming decoupled. If material is prevented from 
reaching the inner disk, the material inside the gap would drain onto the central star on 
a viscous time at ~ 4AU, which with a ~ 10~ 3 , would be ~ 3 x 10 5 yrs. Such a short 
timescale suggests that the probability of observing a system in this stage is very low, but 
still possible. 

Instead of the a ~ 10~ 3 adopted above, Hartmann et al. (1998) estimated a typical 
value of a — 10~ 2 for T Tauri disks based on mm-wave disk mass estimates and assuming 
that these disks expand to ~ 100 AU in ~ 1 Myr. If the TW Hya disk (the region containing 
most of the mass) is not much larger than 200 AU in radius, then a cannot be this large. A 
smaller viscosity could be consistent with the typical T Tauri data if disk masses are a few 
times larger than typically estimated, as suggested by D'Alessio et al. (2000). 

One remaining question is why TW Hya's disk evolution has been so much slower than 
that of other members of the association with presumably similar ages. In part, this must 
be due to the presence of binary (or multiple) stellar companions in several other systems 
(e.g., HD 98800; Soderblom et al. 1998) which can disrupt disks. In addition, we conjecture 
that TW Hya may simply have had a larger initial disk radius (equivalent to a large Ri 
in terms of the similarity solution above), and thus a lower average surface density, than 
other systems. A larger radius would lead to slower viscous evolution for a given a. Also, 
in contrast to timescales for dust coagulation and sedimentation, which depend mostly 
on orbital period (see above), theoretical models of the runaway growth of giant planets 
indicate that the timescale for such growth is a very sensitive function of surface density 
(Pollack et at. 1996). The lower surface density of the disk around TW Hya would imply a 
slower growth of giant planets. This would imply, in turn, that the final clearing of disks is 
due to sweeping-up of material by large bodies. 

6. Conclusions 

The disk of TW Hya seems to be in advanced state of dust evolution. A planet or other 
large perturbing body may have already formed, opening a gap with outer edge around ~ 4 
AU, and although there is still material in the inner disk, it is rapidly dissipating. In the 
outer disk, grains may be reaching the size necessary to start settling towards the midplane. 
This interpretation can consistently explain the large degree of activity of TW Hya, the 
lack of disk emission at near-IR wavelengths and the large fluxes beyond ~ 10 /im, as well 
as the 7 mm images. 
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Evacuated regions in the disks are also known to result from the effects of companion 
stars (GG Tau), and in some cases disk accretion can still occur, especially if the binary 
has an eccentric orbit (DQ Tau). Indeed, HD98800, another member of the TW Hya 
association, shows evidence for an inner disk hole; because this system is quadruple, it 
is quite likely that the inner disk regions have been evacuated by the companion stars. 
Although we cannot rule out such a possibility in TW Hya, the fact that accretion still 
occurs onto the central star, which is not the case in HD 98800, suggests that the companion 
body responsible for the gap is much less massive than a typical companion star, and thus 
much less effective in clearing out the disk. Interferometric imaging should be attempted to 
constrain the mass of any companion object, expected to have a separation of order 2-3 AU 
(0.04 - 0.06 arcsec at 55 pc). 
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Fig. I. — Spectral energy distribution of TW Hya. Observations are from Rucinski & 
Krautter (1983; average of UBVRI measurements) (open circles), Webb et al. (1999) 
(triangles), Sitko et al. (2000) (small open circles), this paper (open square), Jayawardhana 
et al. (1999) (squares), IRAS (de la Reza et al. 1989; Gregorio-Hatem et al. 1987; crosses), 
Weintraub et al. (1989) (hexagons), Wilner (2001) (open circle), and Wilner et al. (2000) 
(filled circles) The solid line is the median SED of classical T Tauri stars in Taurus (Paper 
II). Note the excess flux at mid-infrared wavelengths and the flux deficit below 10 /zm (see 
text). 
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Fig. 2. — SED of TW Hya and disk models with gas and dust well mixed and extending 
to the magnetospheric radius, for a max = 1 mm (solid line), 1 cm (dotted line), and 10 cm 
(long dashed line). Larger maximum grain sizes produce a better match to the sub- mm and 
mm-wave spectrum. The SED of the stellar photosphere is also shown for comparison (short 
dashed line.) 
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Fig. 3. — Model disk adopted for TW Hya. The outer disk, where grains have grown to ~ 
1 cm, has a edge at R ~ 3 — 4 AU, and surrounds the inner optically thin disk, which has 
a vertical optical depth at 10/im t\q ~ 0.05. Gas still exists in the inner disk accreting onto 
the star through a magnetosphere. A minute amount of ~ 1 dust permeates this gas. 
The dotted line is the surface of the outer disk if it extended inward; the resulting SED for 
this model is shown in Figure 4. 
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Fig. 4. — Above: Fit to the TW Hya SED with our composite disk model. Below: Detail 
of the infrared region. The emission from each region is indicated: outer disk (dot-dash), 
edge of outer disk (dashed), inner disk (dotted), star (light solid), total (heavy solid). The 
stellar spectrum is taken from the Allard & Hauschildt (1995) M dwarf library, where a 
model with appropriate effective temperature and gravity (log g =4) could be found. The 
model has been scaled to the observations at K (2.2 /xm). The excess near ~ 4.5/im could 
be CO fundamental emission; the amount of excess (if any) at A < 5/im is very uncertain, 
as it depends critically on the effective temperature (and model) for the stellar photosphere. 
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Fig. 5.— (a) VLA 7 mm images of TW Hya at two resolutions, from Wilner et al. 
(2000). The synthesized beam sizes are 0"82 x 0"61 p.a. 20 (left) and 0"13 x 0"10, p.a. 39 
(right). The contour levels are ±(2, 3, 4, 5, 6) x the rms noise levels of 0.5 and 0.4 mJy/beam. 
Negative contours are dashed, (b ) Simulated VLA 7 mm images of the disk model brightness 
distribution described in the text. The contours and beam sizes are the same as in a. The 
grey scale shows extended emission, most of which remains undetectable with the available 
VLA sensitivity, (c) Difference images obtained by subtracting visibilities derived from the 
model from the VLA 7 mm data. Again, the contours and beam sizes are the same as in a. 
The residual images show no significant signal. 



